arXiv:1502.01983v2 [astro-ph.CO] 18 Nov 2015 


Detecting Primordial B-Modes after Planck 


Paolo Creminelli,*^ Diana Lopez Nacir,'^’'^ Marko Simonovic,*^’^’® 
Gabriele Trevisan^’® and Matias Zaldarriaga® 

^ Abdus Salam International Centre for Theoretical Physics 
Strada Costiera 11, 34151, Trieste, Italy 

^Departamento de Fisica and IFIBA, FCEyN UBA, Facultad de Ciencias Exactas y Naturales, 
Ciudad Universitaria, Pahellon I, 1428 Buenos Aires, Argentina 

'^Institute for Advaneed Study, Einstein Drive, Princeton, NJ 08540, USA 

SISSA, via Bonomea 265, 34136, Trieste, Italy 

® Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 34136, Trieste, Italy 


Abstract 

We update the forecasts for the measurement of the tensor-to-scalar ratio r for various ground- 
based experiments (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloons (EBEX 
10k and Spider) and satellites (CMBPol, COrE and LiteBIRD), taking into account the recent 
Planck data on polarized dust and using a component separation method. The forecasts do not 
change significantly with respect to previous estimates when at least three frequencies are available, 
provided foregrounds can be accurately described by few parameters. We argue that a theoretically 
motivated goal for future experiments is r ~ 2 x 10“^, and that this is achievable if the noise is 
reduced to ~ 1 /xK-arcmin and lensing is reduced to 10% in power. We study the constraints 
experiments will be able to put on the frequency and ^dependence of the tensor signal as a check 
of its primordial origin. Puturistic ground-based and balloon experiments can have good constraints 
on these parameters, even for r ~ 2 x 10“^. For the same value of r, satellites will marginally be able 
to detect the presence of the recombination bump, the most distinctive feature of the primordial 
signal. 


1 Introduction and motivations 

The year 2014 marked the beginning of the R-mode era in cosmology. After the direct detec¬ 
tion of the lensing R-mode signal by Polarbear [1], BICEP2 [2] pushed the constraints on primor¬ 
dial tensor modes using polarization to a level that is competitive with temperature. Given that 
temperature measurements are close to the cosmic-variance limit for the tensor-to-scalar ratio r, 
improvements in the future will practically only come from polarization. Planck [3] measured the 
level of polarized dust emission on the full sky with unprecedented precision (for previous measure¬ 
ments see for example [4, 5]). In this paper we want to look ahead at the future (and futuristic) 
experiments and understand whether the new data on dust polarization substantially change the 
reach expected for the various experiments. 

In looking at the future sensitivity on r, it is important to have in mind some motivated the¬ 
oretical threshold. In the very near future we will explore the region r ~ 0.1, which corresponds 
to simple monomial potentials. If gravitational waves are not detected, is there another motivated 
threshold to reach? One might argue that r ~ 2 x 10“^ is a reasonable goal for future experiments. 
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First of all, it approximately corresponds to the value predicted by potentials that approach asymp¬ 
totically a constant as exp(—(/)/Mp) (Starobinsky model [6], Higgs-inflation [7], etc.). A similar 
number is obtained by looking at the Lyth bound [8] . The excursion of the inflaton during inflation 
is given by 

If one assumes that r monotonically increases going towards the end of inflation, one can conserva¬ 
tively replace r{N) with the one on cosmological scales. The threshold A(/) = Mp then corresponds 
to r = 8N~‘^ ~ 2 X 10“^. A detection of gravitational waves above this level would convincingly 
indicate a trans-Planckian displacement, under the mild assumption that e increases as one moves 
towards the end of inflation. Another way to argue for the same threshold for r is to study the 
consequences of imposing that the scalar tilt is of order 1/N: Ug — 1 = —a/N [9-11]. If this ap¬ 
proximate equality is not an accident, but holds in a parametric window around N = QO, one can 
argue for the existence of a forbidden region in r between 10“^ and 10“^. This second number 
actually depends exponentially on the precise value of the scalar tilt, but 2 x 10“^ corresponds to a 
reasonable lower bound within the present uncertainties on Ug [9]. All these theoretical prejudices 
should be taken with care, but motivate 2 x 10“^ as a relevant figure of merit. Therefore, we will 
look ahead to check which experiments will be able to get to this value of r. 

The outline is rather simple: in Section 2 we explain the method used throughout this paper, 
while in Section 3 we show the result obtained for various experiments. In Section 4 we consider 
more conservative analyses, focussing on possible evidences that the signal is indeed due to tensor 
modes. A similar study was done in Ref. [12] concentrating on the superhorizon R-modes using 

[13]. 

2 Forecasting Method 

2.1 CMB and Noise 

In linear perturbation theory, the coefficients aim of the T-, E- and R-modes of the CMB are 
Gaussian random variables with zero mean and variance 

bill bmm'1 ( 2 ) 


where X,Y = T,E,B. From these, as customary, one defines the “curly” correlators as = 
£{£ + l)C'^^/(27r). Due to parity invariance, only the TT, EE, TE and RR power spectra are 
necessary to characterize the CMB, the others being zero. In our analysis we consider the R-mode 
power spectrum only, so we drop the superscript RR where possible. This is generated by CAMB 

[14] and, since we are solely interested in the forecast for the tensor-to-scalar ratio r, we set all 
cosmological parameters, except r and the optical depth r, to the current best fit values of Planck 

[15] . Although this may look like a rough approximation, r is expected to be only mildly degenerate 
with the other parameters, the biggest degeneracy being the one with r at low multipoles. We 
are going to marginalize over r using a gaussian prior given by Planck analysis [15]. This is a 
conservative approach for satellites since they will have additional information on reionization. On 
the other hand, since large scale polarization measurements are affected by systematics, it is not 
clear how much they will improve the constraints on r. 
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In the presence of white noise due to the instrumentation, the integration over a Gaussian beam 
to go from real space to harmonic space creates a bias A/^ for the estimators of the power spectra. 
At each frequency channel, this is given by [16] 


A/£ 


£(£ + 1 ] 

2tt 


^pix ^pix 




(3) 


where Upix is the noise variance per pixel of size flpix = ©fwhm^ = 0.425 ©fwhm is the 

beam-size variance. Our treatment here is clearly simplistic, since it does not take into account 
systematics. However, these are very experiment-dependent and go beyond the scope of this paper. 


2.2 Foregrounds 

The presence of foregrounds limits our ability in extracting the CMB signal from the data. 
Fortunately, each component scales differently in frequency, and thus it is possible to separate them 
using maps at different frequencies. In this paper we consider two diffuse components: synchrotron 
emission [S) and thermal dust {D). 

The BB power spectra of the Galactic synchrotron emission and thermal dust in antenna 
temperature is given by the following simple parametrization 


As " , 

= {W^fAn , 


W^ = 
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where the parameters are given in Tab. 1. This parametrization fits well the observed power spectra 
[3, 17, 18]. Since for our analysis we are using as a reference the CMB signal, one has to rescale the 
antenna temperature of these components to the thermodynamic temperature of the CMB. The 
conversion is provided by the frequency dependence of the CMB 
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The frequency dependence of synchrotron and dust rescaled to the thermodynamic temperature of 
the CMB then reads 
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Synchrotron emission is the dominant one below 90 GHz^, while dust becomes increasingly 
important for higher frequencies. In our forecasts, for sky coverage of 70% and 20% we have fixed 
the synchrotron amplitude to the value measured by WMAP for the P06 mask [17], while for the 
10% and 1% to the extrapolation of [19] of the WMAP data. The synchrotron spectral index is 
known to steepen as Galactic latitude increases (see e.g. [20]). However we checked that this effect 
can be safely neglected for our purposes (for a detailed study see [21]). 

^The synchrotron and dust power spectra are comparable at roughly 90 GHz for the cleanest 1% of the 
sky. In regions with higher levels of polarized dust emission the transition happens at a lower frequency. 


3 










Parameter 

Synchrotron 

Dust 

^72% [/^K^] 

2.1 X 10-5 

0.169 

^53% [T^K^] 

2.1 X 10"® 

0.065 

^24% [T^K^] 

2.1 X 10-5 

0.019 

^11% [t^k^] 

4.2 X 10-5 

0.013 

7li% [gK^] 

4.2 X 10-5 

0.006 

u [GHz] 

65 

353 

1 

80 

80 

a 

-2.6 

-2.42 

P 

-2.9 

1.59 

T[K] 

- 

19.6 


Table 1: Foreground parameters obtained or extrapolated from [3, 17, 19], as explained in the 
text. refers to the cleanest effective area fsky 

However, polarized dust emission is the leading contaminant for balloon and ground experiments 
probing frequencies higher than 90 GHz. A detailed measurement of the polarized dust has become 
available only recently [3]. For this reason, its impact has been underestimated in some previous 
analyses. For example, the observed value of the dust power spectrum at 353 GHz for 72% of the 
sky is roughly 10 times bigger than the dust-model A used for the forecast of GMBPol [22] . It is now 
clear that there are no regions in the sky for which a measurement of r is achievable without having 
to deal with foregrounds [3, 19]. One of the primary goals of this work is to provide new forecasts 
for the detection and measurement of r with more realistic levels of thermal dust contaminating 
the primordial signal. In this respect, we use the levels of dust presented in [3] for the 353 GHz 
channel of Planck and extrapolate their results when needed. In particular, Planck [3] has recently 
provided the power spectra of dust at 353 GHz for a clean effective area of 72%, 53%, 24%, and 
1%. For the 1%-patch we use as a reference the value of the dust amplitude in the BICEP2 region. 
Even though Planck observed cleaner patches of the same size, the associated errors are too large 
to reliably determine the correct level of dust. 

Some experiments, e.g. Spider, will probe patches of the order of 10% and we need to extrapolate 
the amount of dust on a patch roughly as big as theirs. This can be done in several ways. As a first 
guess, the interpolation of the values provided by Planck as a function of fsky gives = 0.013. 
Another way consists in using the low A^hi region of [23] which covers 11% of the sky. Using 
the relation A'hi = 1-41 x 10^® cm“^ {T 353 ), where T 353 is the optical depth at 353 GHz, we find 
A^hi = 1-35 X 10^^cm“^. Substituting this value in the relation between the amplitude of H-modes 
and A^hi [3], we find that A^^ = 0.013. We tested this procedure against the amplitudes indicated 
by Planck for the 1% of the sky for BIGEP2, and found good agreement. A third way consists in 
using the relation between the amplitude in polarization and the one in intensity [3], Ajy oc 1353. 
This relation for the same low Nhi region gives A^^ = 0.009. We will use = 0.013 as an 
upper bound for the dust levels in the region observed by Spider. 

An additional complication comes from the correlation between synchrotron and dust. It has 
been observed [19] that the correlation among these two components can be as high as 70%. To 
account for this effect in our forecast, for simplicity we will assume that in the power spectrum the 
correlation enters as and set g = 0.5 in our fiducial sky-model, independently of fsky 

and i. 
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2.3 Likelihood and Fisher Analysis 

In the case of an experiment covering the whole sky, one can write the signal d measured at the 
frequency t'i, in harmonic space, as 




(7) 


where W provides the frequency dependence of each component^, the bar indicates that the param¬ 
eters are fixed to their “true” value, and n is the instrumental noise. Assuming that the amplitudes 
are Gaussian (also those of foregrounds), starting from the (and avoiding the use of indices) 


^2 = ^(d -W-a)'^-N-^-{d-W-a) + J- • a, 

£,m 


( 8 ) 


where C is the covariance matrix of the amplitudes of the various components and N is the covari¬ 
ance matrix of the noise, the likelihood can be written as 


CBBid,p) oc I oc e"! (9) 

In order to do a Fisher analysis we are interested in the average of the log-likelihood, 

{log Cbb) = “2 + 1 ) 

I 

( 10 ) 


log Det 


W -W^ +Ni 

W -W^ +Mi 


+ Tr 


W-Cf^-W^+Ni 
W -W^ +Mp ~ 


where the normalization is such that for Cp = Cp, {logCBB) = 0. In the following we will refer to 
the use of this formula as “Component Separation” (CS). 

Given the likelihood as a function of the parameters p, one can define the Fisher matrix 

, ( 11 ) 

p=P 

where p is the set of “true” cosmological parameters. The minimum error on the parameters given 
the data, is set by the Cramer-Rao bound to be 


(log Cbb) 

dpidpj 


4 ^ (12) 

In real experiments, only part of the sky is observed or can be used for cosmology, and thus not 
all modes are available for the analysis. To capture this effect in a simple way one can multiply the 
RHS of Eq. (10) by the fraction of the sky available fsky 

In general the analysis is far more complicated. The likelihood as written in Eq. (10) becomes 
lossy, modes of multipole t will be coupled with their neighbors. In addition, for experiments cover¬ 
ing only a few percent of the sky, the leakage of E- into R-modes will add additional complications 
to the detection of the signal produced by gravitational waves [24] . In our discussion, we will neglect 
these details, and thus the results produced with this method can be seen as optimistic. 

^Since we are considering three components (the index c runs over CMB, Dust and Synchrotron), and 
we are expressing everything in thermodynamic temperature, IT is a 3 x A^channei matrix whose rows are 

{i,wp,wp). 
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In the analysis of Sec. 3, the likelihood in Eq. (10) is a function of 7 parameters^ p = 
{r, Ajj, As, Pd, Psi g,T)y and the error on the tensor-to-scalar ratio is marginalized over the other 
parameters. Notice that at low frequency the effect of changing Td is very similar to a change in 
Pd- As done in [25], we are assuming that the temperature of the dust emitting polarized radia¬ 
tion is the same as the one determined from the intensity maps. In the following we will assume 
gaussian priors for As, Ad, Ps with variance of 50%, 50%, 10% respectively. For Pd we use a 
gaussian prior of 10%, 30% and 50%, for fsky ^ 50%, ~ 5 — 10% and ~ 1% respectively. For the 
optical depth we assume a gaussian prior given by Planck analysis [15]. No priors are assumed for 
the other parameters. 

In Sec. 4 we will instead take a more conservative approach and consider the likelihood in 
Eq. (10) as a function of 11 parameters. Namely, in addition to the 4 previous parameters, we will 
include the ^-dependence of the foregrounds, ud and as, the spectral dependence of the CMB and 
its .^-dependence. In particular we will constrain how much the CMB deviates from what expected 
using simple power laws with parameters acMB for the ^-dependence and Pc mb for the spectral 
dependence. Notice that acMB roughly corresponds to the tensor tilt. 

In all the tables that follow, we use the symbol ” when cTr > r. 


2.4 Likelihood and Fisher Analysis: a phenomenological approach 

The foreground model used in the previous section is of course approximate. For instance the 
foreground parameters may be space-dependent, and the distribution is not exactly gaussian. For 
this reason, and also to compare our results with previous forecasts, we also adopt a second, more 
phenomenological, method for estimating the error on the tensor-to-scalar ratio. It was proposed 
in [26] and already used also in previous forecasts for CMBPol [22] and more recently also in [12]. 
The method assumes that with already known techniques such as [27, 28] it is possible to subtract 
the foregrounds up to a certain level in each single map. The noise introduced by the foreground 
removal is then modeled by accounting the number of possible cross correlations. The sum of 
foreground residual and noise is 




4M 




W. 


Achannel( ^channel 1) 


^tmp 


(13) 


where ap is the fraction of leftover foregrounds in power, lyf includes the conversion to thermo¬ 
dynamic temperature, I'tmp is the reference frequency at which the foreground template has been 
created (e.g. 30 GHz for synchrotron and 353 GHz for dust), and A^channei is the number of fre¬ 
quency channels. This quantity is then treated as an additional source of uncorrelated noise and 
as such is added to the noise bias. The resulting effective noise is given by 

i,j>i ^ ^ 


This method has the advantage of being independent on any specific technique of foreground 
subtraction, but, as already noted in [29] (see also [30]), it has the downside of not being specific 
of any real experiment. By comparing the results of this and the GS method, one can estimate 
what is the level of foreground subtraction obtained by the component separation, and therefore 

^Notice that even in the limit in which there is no additional information coming from the difference in 
^-dependence, an experiment with three frequencies would allow a complete reconstruction of the parameters, 
since its covariance matrix would have six entries. 
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the level at which the foreground modeling needs to be trusted. As already done in [22, 26] and 
recently in [12], in our forecasts we will assume that foregrounds can be subtracted at 1% in power 
level in each map. 

The likelihood in this approach is the one of a single channel with an effective noise bias given 
by Eq. (14), and it reads 


{\og Cbb) 


i 


log ^ 








(15) 


In order to render the forecasts more realistic, we marginalize the error of the tensor-to-scalar 
ratio over the foreground residuals. This can be done by considering the percentage of foreground 
removal as an additional parameter and by multiplying the likelihood in Eq. (15) by a gaussian 
prior describing our ignorance about the exact level at which foregrounds have been removed. We 
will assume that the percentage of foreground removal (7^ is known with relative error of order 1. 


2.5 Delensing 

One of the most important limiting factors in measuring primordial i?-modes is gravitational 
leasing. Since S-modes induced by leasing have the same frequency dependence as primordial ones, 
one cannot proceed in the same way as with foregrounds. The idea then is to use the measurements 
of polarization or temperature on short angular scales to reconstruct the leasing potential and 
remove leasing effects from the i?-mode polarization on the larger angular scales [31-34]. 

Polarization delensing is not always important in improving the constraints on r. In the case 
of a large tensor-to-scalar ratio, the dominant contribution to the signal-to-noise comes from the 
low multipoles. For example, the leased i?-modes power spectrum becomes comparable to the 
power spectrum of the primordial S-modes corresponding to r = O.I around ^ ~ 100. In this case 
delensing cannot improve the errors signihcantly. The other case in which delensing is irrelevant 
is when the noise is larger than the leasing signal. Indeed the power spectrum of leasing B-modes 
for i < 150 is similar to white noise with amplitude 4.4 //K-arcmin, so that only for experiments 
with smaller noise delensing is relevant. If the instrument has a high sensitivity and a good angular 
resolution needed for the reconstruction of the leasing potential, delensing the maps can become 
important and can improve the errors on r even by a factor of a few for the beam sizes and the 
sensitivities of different proposed future experiments [33-35]. For example. Fig. 3 of [34] shows that 
with a beam size of 5 arcmin and sensitivity of I ^uK-arcmin one can improve the constraints on r 
by a factor of 5. Moreover, as the sensitivity and the resolution increase, there are no limits in how 
much of the leasing signal can be subtracted. This improvement is marginal for all the experiments 
considered in this paper, except for the generation-IV experiments (GRD, BAL and ULDB) and for 
the satellites COrE and CMBPol. For all other experiments either the noise level is high enough to 
make delensing irrelevant or the angular resolution is too large to implement delensing. To include 
delensing in those experiments where it is viable, we assume that the power of leasing R-modes is 
reduced to 10% of the original value. The residual power would correspond to a noise equivalent 
power of 1.4 //K-arcmin. It would be useful to explore how the presence of foregrounds impacts the 
ability of delensing, but this goes beyond the scope of our paper. 
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3 Results 

3.1 BICEP2/Keck and Planck 

Let us begin our analysis by checking that our forecasting method is compatible with the recent 
recent joint analysis [36]. In order to do so, we combine the BICEP2/Keck 150 GHz channel with 
the Planck 353 GHz one. With the GS method, for a fiducial value of r = 0.05 using the multipoles 
[20-150] and a gaussian prior on /?£> with variance = 0.11, we obtain an error ar = 0.04, which 
is in good agreement with what was reported in [36]. 

Alternatively, using the phenomenological approach of 2.4, in Fig. 1 we plot the contours rep¬ 
resenting the error on r assuming r = 0 as a function of sensitivity and foreground residuals. As 
can be seen, our ability to constrain r crucially depends on the foreground removal. A reduction 
of foregrounds to 10% in power can lead (in the absence of gravitational waves) to a quite strong 
upper bound r < 0.05 at 2>a. Polarization experiments are already competitive with constraints 
from temperature alone, even with only one frequency and very mild foregrounds removal. Notice 
that with T-modes only cosmic variance prevents to constrain r better than 5 x 10“^ [37].^ Notice 
moreover that in the near future it will be possible to include in the analysis the 95 GHz channel 
of Keck (here we assume a noise of 9 /rK-arcmin). As can be seen in Tab. 2, the inclusion of such a 
frequency would allow to reduce the error on the tensor-to-scalar ratio by a factor of 2 with respect 
to the current constraint. 


r 

BICEP2/Keck -|- Planckssa 

0.1 

3.5 X 10-2 

0.01 

— 

0.001 

— 

0 

2.2 X 10-2 


Table 2: la errors on r for BICEP2/Keck (95 and 150 GHz) and the 353 GHz channel of Planck. 

This error, as calculated from the phenomenological method of Sec. 2.4, corresponds to 
a 30% level of leftover foregrounds in power. 


3.2 Balloon-borne and Ground-based Experiments 

The situation will improve in the next few years since there are several experiments that are 
already collecting data and will have maps in two or more frequencies. In our forecasts we con¬ 
centrate on a few proposed and funded experiments. In particular, for what regards ground-based 
experiments we consider Keck/BICEP3 and the Simons Array, and also AdvACT, CLASS and 
SPT-3G. The specifications used are in Tab. 10. We vary the level of foregrounds according to 
the fraction of the sky targeted by each experiment, as given in Tab. 1. For what concerns the 
available multipoles, we consider the range [30,150]. For AdvACT, CLASS and the Simons Array 
we consider the range [2,150], since they observe a larger fraction of the sky. This range is probably 
larger than what these experiments will actually be able to observe, since to handle atmospheric 
contamination they need to filter the data, losing power at low £’s. In Sec. 4 we will take a more 
conservative perspective. As can be seen in Tab. 3, which summarizes our forecasts, we expect 

^This bound of course is impossible to achieve due to sample variance induced by masking the sky. 







r 

Keck/BICEP3 

Simons Array 

AdvACT 

CLASS 

SPT-3G 


0.1 

1.9 X 10-2 

7.6 X 10-3 

7.3 X 10-3 

6.5 X 10-3 

9.0 X 10-3 

CS 

0.01 

0.001 

9.3 X 10-3 

5.0 X 10-3 

4.5 X 10-3 

3.4 X 10-3 

4.2 X 10-3 


0 

8.1 X 10-3 

3.4 X 10-3 

2.6 X 10-3 

9.0 X 10-^ 

3.7 X 10-3 


0.1 

1.8 X 10-2 

1.1 X 10-2 

9.5 X 10-3 

6.0 X 10-3 

8.1 X 10-3 

EG 1% 

0.01 

0.001 

7.5 X 10-3 

8.1 X 10-3 

6.9 X 10-3 

3.5 X 10-3 

4.1 X 10-3 


0 

6.1 X 10-3 

7.8 X 10-3 

6.6 X 10-3 

3.3 X 10-3 

3.7 X 10-3 


Table 3: la errors on r for future ground-based experiments. 

these experiments to explore values of r of order 10“^. The CS method gives results which are 
roughly comparable to a reduction of foregrounds to 1%. Of course there are sometimes sizable 
differences, since we do not expect all experiments to reduce foregrounds in the same way. Still 1% 
represents a rough estimate of the level at which one should trust the foreground modeling for the 
CS. 



r 

EBEX 10k 

Spider 


0.1 

1.5 X 10-2 

1.8 X 10-2 

CS 

0.01 

7.4 X 10-3 

— 

0.001 

— 

— 


0 

6.4 X 10-3 

1.3 X 10-2 


0.1 

2.2 X 10-2 

2.6 X 10-2 

EG 1% 

0.01 

— 

— 

0.001 

--- 

— 


0 

9.2 X 10-3 

2.1 X 10-2 


Table 4- 1<7 errors on r for future balloon-borne experiments. 

Regarding balloon-borne experiments, we consider EBEX 10k and Spider (which has already 
finished the first flight). As can be seen in Tab. 4, we expect these experiments to explore values 
of r ~ few X 10“^. It is fair to say that the level of dust measured by Planck only slightly degrade 
the previous forecasts and that the goal of these experiments are still within their range. 

Eor all these experiments the error on the dust amplitude will be substantially smaller than the 
present Planck constraints, so that a cross-correlation with Planck will not significantly reduce the 
errors. However, Planck’s data will still be useful to test the spectral dependence of the polarized 
dust emission model. 

Looking a bit further into the future, we also consider an idealized balloon (BAL) and an ultra 
long duration ballon (ULDB) with the same four frequencies of EBEX 10k (150, 220, 280 and 350 
GHz) and beams of 5 arcmin, but leave their sensitivity as a free parameter. For simplicity we 
assume that the sensitivity is equal across the frequencies, even though this may not be the optimal 
choice. The results of our forecasts can be found in Fig. 2 where we estimate the la error for r = 0 
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Figure 1: Contours represent the la error on r for r = 0, marginalized over the foreground 
residuals, using a single map of 1% of the sky at 150 GHz with beam of 30 arcmin, 
as a function of foreground residuals and instrumental sensitivity. This is the case of 
BICEP2/Keck, which has sensitivity of 3-4 ^K-arcmin. 


for BAL covering the few supposedly clean patches found by Planck (~ 5% of the sky), and ULDB 
covering 60% of the sky. As it can be seen, with a noise level ~ 1 /rK-arcmin and lensing removed 
to 10%, which are possible but challenging, one can detect r ~ 2 x 10“^ with high statistical 
significance. Notice however that obtaining a noise level close to 1 /rK-arcmin on 60% of the sky 
with an ULDB seems out of reach with a single 100-days flight (which is the target of this kind of 
experiments). We will discuss about a futuristic ground-based experiment (GRD) in Sec. 4. 


3.3 Satellite Experiments 

Finally, let us present our forecasts for various proposed satellite experiments (see Tab. 11 for the 
specifications). We assume that an effective area of 70% is observed with foregrounds parameters 
given in Tab. 1, and limit the multipole range to [2,300] for GOrE and GMBPol (EPIG-2m) and 
[2,150] for LiteBIRD^. The summary of our estimates can be found in Tab. 5, and it shows that 
with respect to previous forecasts [22] there is only a minor degradation (no more than a factor of 
2) of the ability to detect primordial tensor modes. In particular a detection of r ~ 2 X 10 ^ is still 
achievable for the proposed missions. Notice however that the upper limit for r = 0 below 10“^ are 
very optimistic and degrade significantly when the reionization bump is excluded (we are going to 
comment on this in the next section). 

^The angular resolution of LiteBlRD is not good enough for delensing and there is no advantage in 
considering higher multipoles when lensing is not subtracted. 
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(a) BAL covering 5% of the sky. 


o 

r—H 

X 

b 



(b) ULDB covering 60% of the sky. 


Figure 2: la error on r as a function of the instrumental sensitivity for two hypothetical balloon 
experiments (BAL and ULDB in Tab. 10). These estimates use a patch of 5% of the 
sky for BAL and 60% of the sky for the ULDB. The solid line assumes leasing is not 
subtracted, the dotted line assumes leasing has been reduced to 10% in power. Since 
also delensing is applied, we considered the multipoles [30, 300]. 


4 More conservative analyses 

It is obvious, especially after the case of BICEP2, that a detection of primordial tensor modes 
must convincingly show that the signal is not contaminated by astrophysical foregrounds. If the 
description of foregrounds in terms of few parameters is accurate, we saw that future experiments 
will be able to remove them with very good accuracy. On the other hand, our knowledge of 
astrophysical foregrounds is rather qualitative and it is not clear at what level the model works. 
For example, for r = 2 x 10“^ foregrounds at 150 GHz are larger than the primordial signal by a 
factor of 10 in amplitude at the recombination bump on the cleanest 1% of the sky, and a factor of 
50 in the 70% of the sky. 

There are of course various ways to check that we are observing primordial gravitational waves. 
The primordial signal is homogeneous over the sky and it has Gaussian statistics, contrary to 
what we expect for foregrounds [38]. Other features that are well known about the signal are its 
dependence both in frequency and in i. To study the ability of future experiments to check these 
features, we add to the parameters discussed in the previous section also the possibility of a power- 
law frequency dependence of the CMB signal with vcmb = 150 GHz. Moreover, 

we multiply the tensor mode power spectrum by a power-law ^dependence {(./with 
l^CMB = 80. This roughly corresponds to the tensor tilt, although we are here interested in checking 
the expected approximate scale-invariance and not to assess the possibility to detect the tensor tilt. 
A convincing detection of primordial tensors should constrain both acMB and flcMB to be close to 
zero. This will also give a sense of how close an unmodelled foreground component must be to the 
GMB signal to be undistinguishable from it. Since we want to be more conservative we also add 
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r 

CMBPol 

COrE 

LiteBIRD 


0.1 

1.6 X 10"^ 

1.7 X 10-5 

2.3 X 10-5 

cs 

0.01 

3.0 X 10-^ 

3.5 X 10-^ 

6.7 X 10-^ 

0.001 

1.1 X 10“^ 

1.7 X 10-^ 

3.1 X 10-^ 


0 

2.1 X 10-5 

3.3 X 10-5 

8.9 X 10-5 


0.1 

2.6 X 10-5 

2.6 X 10-5 

3.8 X 10-5 

FG 1% 

0.01 

5.7 X 10-^ 

7.1 X 10-^ 

1.3 X 10-5 

0.001 

3.6 X 10-^ 

5.1 X 10-^ 

1.0 X 10-5 


0 

3.4 X 10-^ 

4.9 X 10-^ 

9.9 X 10-^ 


Table 5: la errors on r for various proposed satellite experiments. For CMBPol and COrE a 
delensing of 10 % has been taken into account. 



r 



^OCMB 

AdvACT 

0.1 

0.001 

1.4 X 10-2 

1.2 X 10-5 

1.8 X 10-5 

CLASS 

0.1 

2.3 X 10-2 

1.9 X 10-5 

2.2 X 10-5 


0.01 




Keck/BICEP3 

0.1 

0.01 

2.8 X 10-2 

1.7 X 10-5 

1.1 X 10-5 

Simons Array 

0.1 

0.01 

1.9 X 10-2 

1.3 X 10-5 

2.2 X 10-5 


0.1 

1.2 X 10-2 

1.3 X 10-2 

2.2 X 10-5 

SPT-3G 

0.01 

7.2 X 10-5 

9.8 X 10-5 

1.1 


0.001 

— 

— 

— 

EBEX 10k 

0.1 

1.6 X 10-2 

4.7 X 10-5 

3.9 X 10-5 


0.01 




Spider 

0.1 

0.01 

3.3 X 10-2 

4.7 X 10-5 

5.4 X 10-5 


Table 6: 1 a errors on r, acMB and jdcMB for future ground-based and balloon-borne experiments. 
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as new parameters the ^-dependence of dust and synchrotron {ao and 05 ) so that the likelihood 
is a function of 10 parameters. 



r 

qrt 

^CXCMB 



0.1 

1.6 X 10-3 

2.2 X 10-2 

6.8 X 10-3 

CMBPol 

0.01 

3.9 X lO-'^ 

5.1 X 10-2 

2.6 X 10-2 


0.001 

2.1 X 10-"^ 

1.4 X 10-3 

1.3 X 10-3 


0.1 

1.7 X 10-3 

2.5 X 10-2 

1.0 X 10-2 

COrE 

0.01 

4.9 X 10-"^ 

6.1 X 10-2 

4.3 X 10-2 


0.001 

2.7 X 10-"^ 

1.7 X 10-3 

2.1 X 10-3 


0.1 

2.7 X 10-3 

3.9 X 10-2 

1.5 X 10-2 

LiteBIRD 

0.01 

1.2 X 10-3 

1.2 X 10-3 

7.3 X 10-2 


0.001 

8.0 X lO-'^ 

3.4 X 10-3 

3.7 X 10-3 


Table 7; 1 a errors on r, acMB o.'^d 13c mb for future satellite experiments. 

The results for ground-based and balloon-borne experiments are reported in Tab. 6 and include 
only values of r for which a significant detection is possible, since only in this case the additional 
parameters acMB and Pcmb are relevant. We see that the next generation of experiments will not 
be able to constrain uqmb and Pcmb, unless r ~ 0 . 1 . 


r 

Simons Array 

AdvACT 

CLASS 

CMBPol 

COrE 

LiteBIRD 

0.1 

1.0 X 10-2 

9.7 X 10-3 

8.3 X 10-3 

2.7 X 10-3 

2.8 X 10-3 

3.7 X 10-3 

pc; 0.01 

0.001 

8.3 X 10-3 

7.6 X 10-3 

6.1 X 10-3 

3.9 X 10-^ 

4.6 X 10-3 

1.0 X 10-3 

— 

— 

— 

1.6 X 10-3 

2.4 X 10-3 

7.5 X 10-3 

0 

8.1 X 10-3 

7.4 X 10-3 

5.9 X 10-3 

1.4 X 10-3 

2.1 X 10-3 

7.2 X 10-3 


Table 8: la errors on r for big-patch experiments, assuming I > 30. 

Our results for satellite experiments are given in Tab. 7: we find that the inclusion of additional 
parameters does not significantly degrade the errors on r (at most by a factor of 2) . Even for r = 0.1 
the check of the tensor consistency relation, which would give acMB — 10 “^, looks impossible. 

Another point of concern about foregrounds is the possibility of detecting the reionization bump. 
This of course is only relevant for experiments looking at a large portion of the sky. At this stage 
our knowledge of polarized foregrounds on large scales is very limited and it is not clear whether 
the reionization bump will be accessible once foregrounds are included. Moreover, ground-based 
experiments will also be limited by atmospheric contaminants. While in the previous sections we 
extended the analysis to low multipoles for experiments with large fsky, in Tab. 8 we consider a more 
conservative analysis where only the multipoles £ > 30 are considered. To do so we consider the 
likelihood a function of six parameters, as in the previous section. While the change is moderate for 
large values of r, since the low multipoles do not help much with the statistics, the effect is relevant 
for small values of r and becomes dramatic for r = 0 ; the upper limit is degraded by a factor of 
10 (this is compatible with the results of [34]). Notice that the amplitude of the reionization bump 
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depends strongly on r, so that when the measurement of tensor modes relies on the large scales, 
the error on r is significantly affected by the uncertainty on r. 

Let us now comment on our fiducial threshold r = 2 x 10“^. Ground-based experiments of 
the so-called stage IV are expected to achieve a sensitivity of the order 1 /rK-arcmin with 0(10^) 
detectors over 5 years. For this value of r, in Fig. 3 we show the error on r, acMB and 13c mb for 
a hypothetical ground-based experiment (GRD) as a function of fsky for two different sensitivities. 
The detection of r = 2 x 10“^ can be achieved at more than 3(T if the maps are delensed to 10% 
and roughly 20% of the sky is observed. In this case the constraints on acMB and 13c mb are small 
enough to allow a clear distinction from our modeled foregrounds. For satellite experiments, from 
our results shown in Tab. 7 and 8, we see that r = 2 x 10“^ is still detectable with large significance, 
even when it is not possible to detect the reionization bump. From Tab. 7 we also see that the 
error on 13c mb is small enough to allow for a clear distinction of 13c mb from /3d (or 13s)- 



0 0.1 0.2 0.3 0.4 0.5 


/sky 


Figure 3: Top panel: Error on r = 2 x 10“^ as a function of fgky for a hypothetical ground 
experiment (GRD in Tab. 10). Solid lines assumes leasing is not subtracted, dotted 
lines assumes leasing has been reduced to 10%. Since also delensing is applied, we 
considered the multipoles [30, 300]. 

Middle panel: Error on acMB- 

Bottom panel: Error on f3cMB- Notice that only for fsky ^ 30% the constraints on 
the spectral dependence allows to reject dust at 2a. 
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The skeptical reader may be worried about the possibility of detecting gravitational waves 
buried under a foreground signal: we can model foregrounds, but how can we be sure that what is 
left in map is due to tensor modes and not some additional “evil dust” component we are unaware 
of? It is fair to say that a robust detection of primordial tensor modes requires a detection of the 
recombination bump. This feature, like a resonance in particle physics, should be robust against 
foregrounds which are not expected to peak al £ ~ 80. To assess the ability of future experiments 
to measure the bump and distinguish it from a featureless power-law dependence, we compare the 
analysis in Sec. 3 (extended to include ucmb): with a model which does not include the tensor 
transfer function, so that the spectrum is just a power law in L By treating the amplitude of the 
bump as a continues parameter, one can use Wilks’s theorem.® In Tab. 9 we report for some of the 
experiments the minimum value of r for which a 3 cj evidence of the bump (compared to a featureless 
power law .^-dependence) is possible. The rule of thumb is that a lOci measurement of r gives a 3(T 
evidence of the recombination bump. For r = 2 x 10“^ it will be challenging to obtain evidence 
of the recombination bump even with satellite experiments. Indeed, for such a small value of r 
the primordial power spectrum at the recombination peak is comparable to the leasing B-modes 
reduced to 10% in power or 1.4 //K-arcmin. 


Simons Array 

AdvACT 

CLASS 

GRD 

CMBPol 

Train 0.08 

0.055 

0.095 

0.005 

0.003 


Table 9: Minimum value of r for which a 3a detection of the recombination bump is possible. For 
GRD we choose a noise of 1 /rK-arcmin, 10% delensing and 20% of the sky. Since we 
are interested in the recombination bump, the analysis is restricted to i > 30. 


5 Conclusions 

In this paper we updated the forecasts for various future B-mode experiments taking into 
account Planck data on foregrounds. For experiments with at least three frequencies, the forecasts 
on r do not change significantly with respect to previous estimates, provided a simple modeling of 
foregrounds in terms of few parameters works at the required accuracy. 

In particular we focussed on the theoretically motivated target of r = 2 x 10“^. This is 
achievable both with balloon-borne and ground based experiments if the noise can be reduced to 
~ 1 /rK-arcmin and lensing B-modes are reduced to 10%. The ground-based experiments covering 
> 30% of the sky should also have the statistical significance to check that the gravitational wave 
signal has a frequency dependence compatible with the one of the CMB and very different from 
the known foregrounds. Even for satellite experiments observing the recombination bump, which 
would likely be a convincing evidence that the signal is indeed due to primordial tensor modes, will 
be challenging for r = 2 x 10“^. 

®We compute the minimum value of r for which | (1 — CDF (2(log£bump) — 2(log Tnobump))) < 0.003, 
which corresponds to a 3a confidence level, where CDF is the cumulative distribution function for the yf 
distribution with one degree of freedom. 


15 






Acknowledgements 

It is a pleasure to thank C. Baccigalupi, D. Gaggero and F. Paci for useful discussions. We 
would also like to thank J. Errard and G. Fabbian for pointing out the new specifications of the 
Simons Array. We also would like to thank the unknown referee for the careful review of our paper. 
M.S. gratefully acknowledges support from the Institute for Advanced Study. M.Z. is supported in 
part by NSF Grants No. PHY-1213563 and No. AST-1409709. 


16 



A 


Instrumental Specifications 


Experiment 

fsky [%] 

ly [GHz] 

dpWHM n 

dP [pK'] 



90 

2.2 

11 

AdvACT 

50 

150 

1.3 

9.8 



230 

0.9 

35 



38 

90 

39 



93 

40 

13 

CLASS 

70 






148 

24 

15 



217 

18 

43 



150 

6.6 

5.5 



220 

4.7 

11 

EBEX lUk 

2 






280 

3.9 

25 



350 

3.3 

52 



95 

30 

3.0 

Keck/BICEP3 

1 

150 

30 

3.0 



220 

30 

10 



95 

5.2 

13.9 

Simons Array 

65 

150 

3.5 

11.4 



220 

2.7 

30.1 



94 

49 

17.8 

Spider 

7.5 

150 

30 

13.6 



280 

17 

52.6 



95 

1 

6.0 

SPT-3G 

6 

150 

1 

3.5 



220 

1 

6.0 

BAL 

5 

150, 220, 280, 350 

5 

[1,5] 

ULDB 

60 

150, 220, 280, 350 

5 

[1,10] 

GRD 

[1,50] 

100, 150, 220 

5 

1, 3 


Table 10: 


Specifications for balloon-borne and ground-based experiments used in our forecasts: 
[39-4-6]. The sensitivity 5P = (JpixOpwHM is for the Stokes Q and U. 
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Experiment 

u [GHz] 

dpWHM n 

SP [fiK'] 


30 

26 

19.2 


45 

17 

8.3 


70 

11 

4.2 

CMBPol (EPIC-2m) 

100 

8 

3.2 


150 

5 

3.1 


220 

3.5 

4.8 


340 

2.3 

21.6 


45 

23 

9.1 


75 

14 

4.7 


105 

10 

4.6 


135 

7.8 

4.6 


165 

6.4 

4.6 

COrE 

195 

5.4 

4.5 


225 

4.7 

4.6 


255 

4.1 

10.5 


285 

3.7 

17.4 


315 

3.3 

46.6 


375 

2.8 

119 


60 

32 

10.3 


78 

58 

6.5 

LiteBIRD 

100 

45 

4.7 

140 

32 

3.7 


195 

24 

3.1 


280 

16 

3.8 


Table 11: Specifications for satellite experiments used in our forecasts: [22, fl, f8]. The sensitivity 
6P = (JpixOpwHM is for the Stokes Q and U. All experiments target approximately 70% 
of the sky. 
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